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In this paper polarized valon distribution is derived from unpolarized valon distribution. In driving polarized valon distribution some 
unknown parameters exist which must be determined by fitting to experimental data. Here we have used Bernstein polynomial method 
to fit QCD predictions for the moments of structure function, to suitably the constructed appropriate average quantities of the El 43 
and SMC experimental data. After calculating polarized valon distributions and all parton distributions in a valon, polarized parton 
density in a proton are available. The results are used to evaluate the spin components of proton. It turns out that the results of polarized 
structure function are in good agreement with all available experimental data on g^" of proton. 



1 INTRODUCTION 

In valon model, the hadron is envisaged as a bound state 
of valence quark cluster called valon. For example, the 
bound state of proton consists of two U and one D valons. 
These valons thus bear the quantum numbers of respective 
valence quarks. Hwa [ 1 1 found the evidence of valons 
in the deep inelastic neutrino scattering data , suggested 
their existence and applied it to a variety of phenomena. 
Hwa and Yang published two papers [ 1 1 and improved 
the idea of valon model and extracted new results for the 
valon distributions. Recently we published a paper [ |2| 
in which we extracted unpolarized constituent quark and 
hadronic structure function in Next-to-Leading order. We 
were interested here to extend valon model to polarized 
ones which has not yet been done. 



2 Unpolarized and Polarized Valon Distri- 
butions 

Valence and its associated sea quarks plus gluons in the 
dressing processes of QCD is defined as valon. In a bound 
state problem these processes are virtual and a good ap- 
proximation for the problem is to consider a valon as an 
integral unit whose internal structure cannot be resolved. 
The proton, for example, has three valons which, on the 
one hand, interacts with each other in a way that is charac- 
terized by the valon wave function and which on the other 



hand responds independently in an inclusive hard collision 
with a dependence that can be calculated in QCD at 
high Q^. This picture suggests that the structure function of 
a hadron involves a convolution of two distributions: valon 
distribution in proton and parton distributions in a valon. 
In an unpolarized situation we may write: 



y 



(1) 



the summation is over the three valons. Here F'^ix, Q^) 
is a proton structure function and is the corresponding 
structure function of a v valon and Gv/p{y) indicates the 
probability for the v valon to have momentum fraction y 
in the proton. In Ref. [1] the unpolarized valon distribu- 
tion in proton with a new set of parameters for Gu/p(y) and 
Go/piy) has been recalculated. Now we can construct the 
polarized valon distributions in proton which are based on 
the definitions of unpolarized valon distributions. The po- 
larized valon is defined to be a dressed polarized valence 
quark in QCD with the cloud of polarized gluons and sea 
quarks which can be resolved by high probes. In [ 3 1 
we found the following constrain for the polarized parton 
distributions, \6f{x, Q^)\, and unpolarized one, f{x, Q^), at 
low value of which by positivity requirements implying 



mx,Q^)\<f(x,Q^) 



where / 
rules 



u,u,d,d, s, s,g, and furthermore by the sum 
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Au + Au-Ad + Ad^Ai^ 1.2573 + 0.0028 , 
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Au + Au + Ad + Ad- 2{As + As) = Ag = 0.579 + 0.025 



3 Moment otg'! polarized structure function 



here Af is the first (n - I) moment of Sf, which is defined 
by 



Jo 



Af(Q^)^ I dx6f{x,Q^) 



(2) 



For determination polarized parton distributions, the fun- 
damental analysis is to relate the polarized input density to 
unpolarized ones using some intuitive theoretical argument 
as guide lines. So we can introduce the following equations 
to relate the Non-singlet(NS) and Singlet(S) polarized val- 
ons to unpolarized valon distributions 

dGf/(y) = 26Gu/piy) + SGD/p(y) 

= 26Fu(y) X Gu/p(y) + 6Fo(y) x Go/piy) , (3) 
6Gl(y) = 2f(y) 6Gu/p(y) + f(y) SGo/piy) 

= 2f(y) 6Fu(y)xGu/p(y)+f(y) 6FD(y)xGD/p(y) ,(4-) 

where 6Fj(y) and f(y) are defined by: 

6Fj(y) = Njy^Hl - yf^l + yjy + rjjy'^') , (5) 

f(y) = Ky'-' +Ay+ fxy'' + vy^ + py^' + Ty\ (6) 

the subscript j refer to U, D, and Gu/p(y), Go/piy) has been 
defined in Ref.[l]. The factor 2 in Eqs.(3,4) backs to ex- 
istence of 2-U type valons. By using experimental data 
for g'^ [j4j and using Bernstein polynomials which we will 
explain them in Sec. 4, we did a fitting, and could get pa- 
rameters of Eqs. (5,6) which are defined by unpolarized 
valon distributions U and D in Eqs. (3,4). According to 
the Eq. (2), we can write as an example, the first moment 
of polarized m- valence quark in valon model as following: 

Am,, - I dxSuy - 2 I dy6Gu/p(y), 
Jo Jo 

therefore we can construct for polarized valon distribution 
the following sum rules 

2 I dySGuipiy) - I dydGo/piy) = A3 , 
Jo Jo 

2 I dySGu/p(y) + I dydGo/piy) = Ag . 
Jo Jo 

Here we have not considered the SU(3) symmetry breaking 
where with totally flavor-symmetric we have for light sea 
density 6u - 6d - 6s. So what we will get for A3 and 
Ag values will have a little bit different with respect to the 
quoted values in the above. 



Let us define the Mellin moments of any structure function 
h{x, Q^) as following: 



M{n,Q^)= x"-^hix,Q^)dx . 
Jo 



(7) 



Correspondingly in n-moment space we indicate the mo- 
ments of polarized valon distributions for NS and S sector, 
AGf(n) and AG^(«) as: 



AGf'^(«)= ff-'6G''/''(y)dy. 
Jo 



(8) 



For the moments of polarized singlet and non singlet dis- 
tributions we shall use, the leading order solutions of the 
renormalization group equation in QCD. They can be ex- 
pressed entirely in terms of the evolution parameter s 



In Q^/A^ 
lne^/A2 



5 = In 



(9) 



where 2o, A are scale parameters and we fixed it by = 
IGeV^ and used A = 0.203GeV. It is completely known 
the moments of the unpolarized Singlet and Non-singlet 
and for polarized ones we have 



AM^^in, Q^) = exp(-6d'"'''s) 



<0)n ^ 



(10) 



AM^ (n, g-^) = -(1 + Sp) exp{-6d+s) 



+ -(I - 6p)exp(-dd_s) 



(11) 



where 5p and other associated parameters are as follows: 



6p = 



6/^l"-6df^"+4f6df^" 



A — 6d+ — 6d- 



6d^ = ^-i6d^°l"+6df^"±A), 
\2n 



(12) 



/ is the number of active flavors. The anomalous dimen- 
sions which are simply the n-th moment of polarized 
LO splitting function are given by [ 5 1 . The moments of all 
polarized u and li- valence quark in a proton are: 



AM„„(«, eO + AMrf^,(«, = AM'^^(«, e2)x AGf («)(13) 
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The moment of polarized singlet distribution in a proton 
(2) is as follows: 



n-k 



In Eq.(14) S symbol indicates Yjq=ii,d.s(l + '^^is situa- 
tion we can compute the moments of polarized sea quarks 
as: 



2/ 



(15) 



In moment space we can get the polarized moment of pro- 
ton in terms of their polarized constituent quarks as below 



M(n, e^) = i ^ (e^)(AMg + AM^) . 



(16) 



q=u,d,s 



By inserting Eqs. (13-15) in above equation, 16 unknown 
parameters will be appeared which should be determined 
by a fitting processes. 

4 QCD fits to average of moments using 
Bernstein polynomials 

Because for a given value of Q^, only a limited number of 
experimental points, covering a partial range of values x, 
are available, one can not use the moments equation like 
the present one. A method device to deal to this situation is 
that to take averages of structure functions with Bernstein 
polynomials. 

We define these polynomials as 



Pn,k(x) 



T(n + 2) 



T{k + \)T{n + k + \) 



(17) 



(-1)' 



(n-k-l) 



■M((^ + 0+ 1,2') . 



(19) 



^^^^ where moments M are given by 



M{{k + l) + \,Q^) 



Jo 



(jt+/+i)-i 



gi(x,Q )dx , 



(20) 



The integral (18) represents an average of the function 
gi(x) in the region x„_jt - jAxn^k<x<x„^t + jAx„k where 
Xn,k is the average of x which is very near to the maximum 
of Pn^kix), and Ax„ i is the spread of x„^k- The key point 
is that values of g\ outside this interval contribute little to 
the integral (18), as Pn,k(x) decreases to zero very quickly. 
So, by suitably choosing n, k, we manage to adjust the re- 
gion where the average is peaked to that in which we have 
experimental data. This means that we can use only 41 
averages g„^k like : 



82. 



''^'\Q\g'^\Q\ 



Other restriction which we assume here, is that to ignore 
the effects of moments with high order n which are not very 
effective. To obtain these experimental averages from the 
E143 and SMC data for xgi [ 4|, we fit xgx{x, Q^) for each 
bin in separately, to the convenient phenomenological 
expression 



xg[ 



phen) 



■ Jlx'^il 



xf 



(21) 



this form ensures zero values for at x = 0, and 
X - I. Using Eq.(21) with the fitted values of J[,S,C, 
one can then compute g^"'^*(6') using Eq.(18), in terms of 
Gamma functions. Using Eq.(19) the Bernstein averages 
gn,k{Q^) can be written in terms of odd and even moments 
M{n, Q\ 



52,i(e') = 6(M(2,e2)_A4(3,g2)) 

^2,2(e') = 3(M(3,e2)) , 



Thus we can compare theoretical predictions with exper- 
imental results for the Bernstein averages, which are de- 
fined by [0 



?,a(2')= I dxp„k(x)gi(x,Q^) , 
Jo 



(18) 



In Eq. (17), p„,k(x) are normalized to unity, dxp„ k{x) — 
1. Using the binomial expansion in Eq.(17), it follows that 
the averages of g\ with p„^k{x) as weight functions, can be 
obtained in terms of odd and even moments. 



?n,k 



(n-k)\T(n + 2) 
T{k+ \)Tin-k+l) 



gmo(Q^) = 4004.000001 (M(11, e') - M(14, Q^)) , 
-12012 (M(12, e^) - M(13, e^)) . (22) 

We shall use the result of Eq.(21) for the QCD prediction 
of M{n, Q^). Thus there are 16 parameters to be simultane- 
ously fitted to the experimental g„_(t(2^) averages. Defining 
a global for all the experimental data points, we found 
an acceptable fit with minimum ;^f^/d.o.f. = 0.03262/123. 
Since polarized valon distribution are now determined with 
16 known parameters, the Eqs. (5,6) are now completely 
determined for polarized NS and S valon distributions. In 
Fig.(l) we plotted y x SG^,^ and y x 5Gp as a function of y. 
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Figure 1. Non-singlet and Singlet polarized valon distributions 
as a function of y. 

5 Polarized parton distributions and 
structure function 

Now we want to compute the polarized structure func- 
tion of proton. Since we calculated before polarized valon 
distribution in a proton, by having the polarized structure 
function in a valon, it is possible to extract polarized parton 
structure in a proton. To obtain the z dependence of parton 
distributions in practical purposes from the n-dependent 
exact analytical solutions in Mellin-moment space, one 
has to perform a numerical integral in order to invert the 
Mellin-transformation. Consequently we can get the fol- 
lowing expressions for polarized parton distributions in a 
proton: 

du„(x,Q') ^ 2f 5f'{-,Q^)x6Guipiy)—, 

Jx y y 
Jx y y 

f 1 X 

6Ux,Q') = 6f(-,Q^)xf(y)(2SGu/p(y) 

Jx y 

+6GD/p(y))— , (23) 

y 

In leading order approximation in QCD, according to the 
quark model, g'^ can be written as a linear combination of 
6q and Sq, 

where are the electric charges of the (light) quark-flavors 
q - u,d, s. We are now in a position to present the results 
for the proton polarized structure function, g'^ We pre- 
sented in Fig. (2) the results of g^ix, Q^), for - 2 GeV^ 
and compared it with experimental data [4] . 




X 

Figure 2. The results of g'^ix, Q~), at 2^ = 2 GeV'. Our model 
is indicated by solid line. Other theoretical predictions are from 
Ref.[3,7]. Experimental data points are from Ref.[3] 
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